Thermoelectric effects of an Aharonov-Bohm interferometer with an embedded 

quantum dot in the Kondo regime 
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Thermoelectric effects are studied in an Aharonov-Bohm (AB) interferometer with an embed- 
ded quantum dot in the Kondo regime. The AB flux-dependent transmission probability has an 
asymmetrical shape arising from the Fano interference between the direct tunneling path and the 
Kondo-resonant tunneling path through a quantum dot. The sign and magnitude of thermopower 
can be modulated by the AB flux and the direct tunneling amplitude. In addition, the thermopower 
is anomalously enhanced by the Kondo correlation in the quantum dot near the Kondo temperature 
(Tk)- The Kondo correlation in the quantum dot also leads to crossover behavior in diagonal trans- 
port coefficients as a function of temperature. The amplitude of an AB oscillation in electric and 
thermal conductances is small at temperatures far above Tk, but becomes enhanced as the system 
is cooled below Tk- The AB oscillation is strong in the thermopower and the Lorenz number within 
the crossover region near the Kondo temperature. 

PACS numbers: 72.10.Fk, 72.15.Qm, 03.65.Ta, 73.63.Kv 



I. INTRODUCTION 

The phase-coherence of an electron's wave function in 
mesoscopic systems enables one to observe a wide vari- 
ety of interesting phenomena in quantum physics. One 
of them is the phase shift experienced by electrons due 
to scattering centers. Although the phase cannot be 
measured directly in bulk systems, mesoscopic systems 
provide an opportunity to attain such information. Re- 
cently the transmission phase shift experienced by elec- 
trons passing through a quantum dot was measured us- 
ing the Aharonov-Bohm (AB) interferometersii2i2ii. In 
addition to the direct observation of the phase, phase- 
coherence leads to such phenomena as the magnetic and 
electrostatic AB effects, Fano interference, persistent cur- 
rents, universal conductance fluctuations, etc. 

Quantum dots provide a unique opportunity to study 
nonequilibrium many-body effects due to the strong 
Coulomb interaction and quantization of energy levels, 
spin, and charge. A Fano resonance was observed 5 - in the 
conductance in transport through a quantum dot which 
is connected to two leads. The opening of direct tunnel- 
ing between two leads was believed to be the main reason 
for the observation of the Fano interference pattern. The 
interference between two current paths, resonant and di- 
rect tunneling, leads to the asymmetrical Fano resonance 
in the differential conductance. However, the nature of 
the direct tunneling path is not clear and the direct tun- 
neling cannot be controlled in this work 5 . Inserting a 
quantum dot in an AB interferometer, the Fano shape 
in the IV curves was also observed^. Free control of the 
direct tunneling probability is possible using a gate volt- 
age. 

In addition, the asymmetrical shape of a Fano reso- 
nance has been observed in other nanoscopic systems. 
When magnetic adatoms are deposited on a metallic sur- 
face, the local electronic density of states (DOS) on a 



metallic surface is modified due to the Kondo effect oc- 
curring at the site of a magnetic adatom. Using a scan- 
ning tunneling microscope (STM)L&, the differential con- 
ductance between the STM tip and the metallic surface 
was observed to be of asymmetrical Fano resonance type 
close to the magnetic adatom at low temperature. Con- 
fining electrons to an ellipse 9 in a quantum corral, the 
phase coherence of an electron's wave functions was de- 
tected by imaging the metallic surface with STM spec- 
troscopy. When a magnetic adatom was positioned at 
one of the two foci of the ellipse, the STM experiments 
showed the coherently reproduced image of the magnetic 
atom at the other empty focus point. A Fano interfer- 
ence pattern was also observed in the carbon nanotubesi^ 
when Co ions were deposited on the surface of carbon 
nanotubes. The ability to observe the Fano resonance 
shape in all these experiments depends on the phase co- 
herence of an electron's wave functions propagating along 
a metallic surface. 

In this paper, we consider a phase-coherent Aharonov- 
Bohm interferometer with an embedded quantum dot 
(see Fig. ^| and study theoretically the thermoelectric 
effects of this systemii when the embedded quantum 
dot lies in the Kondo regime. The quantum dot in 
the Kondo regimsiSiiiiii 5 . will be called a Kondo dot. 
Though the issue of the transmission phase in this sys- 
tem is a very interesting subject, we focus on the effects 
of phase-coherence on the transport coefficients. The 
phase-coherence enables us to study the effects of the 
Aharonov-Bohm oscillations as well as the Fano interfer- 
ence on the thermoelectric properties of the model sys- 
tem. We use the nonequilibrium Keldysh Green's func- 
tion method and the noncrossing approximation for our 
model study. 

There are two possible paths for the flow of electrons 
in our AB interferometer. The direct tunneling ampli- 
tude between two leads is a temperature-independent 
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FIG. 1: Schematic display of Aharonov-Bohm(AB) interfer- 
ometer. Electrons can flow from one lead to the other through 
two paths, direct tunneling (Tlr) and a resonant tunneling 
via a quantum dot (VdL, V dR ). The magnetic flux $ is thread- 
ing through the AB ring formed of two current paths. The 
AB phase cj> = 27r«l>/ ( I>o is embedded in the tunneling matrices 
by the overall phase of V dL T LR V Rd = \V dL T LR V Rd \ e^. 



constant, while the tunneling through the Kondo dot is 
strongly temperature dependent due to the Kondo ef- 
fect. Different physics is expected above and below the 
Kondo temperature (Tr-), due to the crossover behavior 
in the Kondo dot as a function of temperature. Since 
the flow of electrons through a quantum dot is blocked 
by the strong repulsive Coulomb interaction (Coulomb 
blockade) at high temperatures above Tr-, practically no 
Fano interference arises and the AB oscillations are very 
weak in transport coefficients. With lowering tempera- 
ture below Tk, the Kondo resonance peak in a quantum 
dot develops close to the Fermi level. This opening of 
an additional resonant current path leads to the Fano 
interference with the direct tunneling path. Fano inter- 
ference and AB flux modify the shape of the transmission 
probability through the system near the Fermi level. 

Out of all the transport coefficients, thermopower is 
one of the most sensitive to the shape or the particle- 
hole asymmetry of the transmission probability and is an 
appropriate experimental probe to investigate the Fano 
interference and the AB effect. We find that the sign and 
magnitude of the thermopower can be modulated by con- 
trolling the Aharonov-Bohm phase with magnetic fields 
and the tunneling matrices with varying gate voltages. 
In addition, we find that the Kondo correlation enhances 
the amplitude of AB oscillations in diagonal transport co- 
efficients at low T compared to T&-. Since two diagonal 
transport coefficients, the electric and thermal conduc- 
tances, are influenced by Fano interference and AB flux 
in the same way, their ratio, the Lorenz number, is insen- 
sitive to AB flux near the zero temperature and is fixed at 
the Sommerfeld value. The amplitude of AB oscillations 
in the thermopower and Lorenz number is strong near 
the Kondo temperature. A short paper which presents 
some of these results appeared elsewhere^. 

This paper is organized as follows. In Sec. [HJ the 
model Hamiltonian is introduced and the formulation of 
electric and heat currents is summarized leaving the de- 
tails to Appendix and [5] I n Sec. IIIII the modified 



noncrossing approximation (NCA) is briefly introduced 
for our model system. Numerical results from solving the 
NCA equations self-consistently are presented in Sec. llVl 
and a conclusion is included in Sec. [3 



II. 



MODEL AND FORMALISM 



In this section, we formulate the electric and heat cur- 
rents in the Aharonov-Bohm interferometer which con- 
tains a quantum dot in one arm. The model system is 
shown in Fig. ^ Using the Keldysh's nonequilibrium 
Green's function methocUiiS, the currents are derived in 
a Landauer-Buttiker form. The phase-coherence of the 
model system enables us to study the Kondo effects of 
a quantum dot, the Fano interference between the direct 
path and the Kondo-resonant path, and the Aharonov- 
Bohm effects by applying magnetic fields. 

When the quantum dot lies in the Kondo regime or 
when there is an odd number of electrons within a dot 
and the spin of a highest-lying electron is unpaired, the 
quantum dot can be effectively considered as a magnetic 
impurity with spin i . The highest-lying electron can then 
be described by the Anderson impurity model treating 
other electrons as a background. The model Hamiltonian 
may be written as H = Hq + Hi where 



Ho = T Ye rcK c t 

p=L,R l a 

+e d d l d a + Un ini , 

o=T,l 

yE [ TLRC lk a c Rk- a + H - c 

kk' a. 



(la) 



Hi 



p=L,R kc 



(lb) 



Hq describes the isolated system of a dot and leads and 
the tunneling between a dot and two leads is described by 
the tunneling Hamiltonian Hi . Here c p g Q is the electron 
annihilation operator of spin direction a in the p = L 
(left) p = R (right) lead. d a is the annihilation opera- 
tor of the highest-lying electron with unpaired spin in a 
quantum dot. Tlr is the direct tunneling matrix from 
the right lead to the left one and Vd p is the hopping ampli- 
tude from the left or right electrode to a quantum dot(d). 
In our model study, the wave vector dependence of these 
tunneling matrices is neglected. 

Since electrons can flow from left to right in two dif- 
ferent paths (direct tunneling and resonant tunneling 
via a quantum dot), quantum interference is expected 
when the phase coherence of electrons is retained. Ap- 
plying a magnetic field to the system, we can study the 
Aharonov-Bohm(AB) effect on the transport. The AB 
phase (j> — 27T$/$ is incorporated into the tunneling ma- 
trices in a manner that VcilTlrYrl = \V d L TlrVrl |e 10 . 
<f> is the magnetic flux passing through the system as 
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shown in Fig.^ an d &o = hc/e is the flux quantum. Ap- 
plied magnetic fields are assumed to be not strong enough 
to split the doubly degenerate spin states of a quantum 
dot. 

The electric current operator can be defined as a 
change in the number of electrons per unit time in the 
left electrode. 
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Here Nl = X)fc Q c Lk ~ c Lka * s ^ ne num ber operator in the 
left lead. The heat current operator is also defined as a 
change in the thermal energy per unit time in the left 
lead. 



Qi 



= ~[H L ,H] 
in 



(3) 



H L is the Hamiltonian of the left lead without the chemi- 
cal potential shift due to the source-drain voltage. Using 
the current conservation in a steady state, both currents 
can be written in a Landauer-Biittiker form in terms of 
the Green's function of a dot. 



Il 
Ql 



de 



T(e) [j^e) -/ R (e)]. (4) 



Here f p (e) — f(e — jj, p ,T p ) is the Fermi-Dirac thermal 
distribution function in the lead p = L, R when each lead 
is in thermal equilibrium at temperature T p . fi p = —eV p 
is the chemical potential shift due to the applied source- 
drain bias voltage. T(e) is the transmission probability 
spectral function through an AB ring and is related to 
the Green's function G r d of a dot by the equation, 



T(e) = T (e)-Im[Z LR G r d }. 



(5) 



The first term To(e) comes from the direct tunneling and 
is given by the expression. 



T (e) 



4jN L (e)N R (e) 
\Dr(e)\* • 



(6) 



Here 7 = ^NlNrITlrI 2 is the dimensionless measure 
of direct tunneling of electrons between two leads. Nl^r 
is the density of states (DOS) at the Fermi level of the 
left, right leads, respectively. The overlined DOS is nor- 
malized such that its value is unity at the Fermi energy. 
D r (e) — 1 — 70i(e)<7jj(e) an d the definitions of reduced 
Green's functions gl' a (e) of each lead (p — L, R) can be 
found in Appendix ^ The resonant tunneling through 
the quantum dot and the interference effect are all in- 
cluded in the second term of T(e). 



jlr 



4Z r L Z R 



rL(e) + r fl ( e ) 



(7) 



In the absence of the direct tunneling, Z^ R is reduced 
to the familiar expression, AT R / (T l + T R ). The pa- 
rameters of Z r L R and r^^e) are tabulated in the ta- 
ble [J Detailed derivations and definitions of new nota- 
tions can be found in Appendix [S] The line-broadening 
parameter T p — irN p \Vd p \ 2 (p = L, R) measures the hop- 
ping rate of electrons between the quantum dot and the 
leads. The complex number z — k 2 N lN R VdLTLRYRd = 
tJ^TlXr e 1 ^ contains the effect of applied magnetic fields 
through the Aharonov-Bohm phase <f> = 27r$/$ - This 
combination of tunneling matrices VdLTL R V R d has the 
meaning that electrons hop off the quantum dot into 
the right lead, tunnel from the right lead to the left one 
via the direct tunneling, and hop back on the quantum 
dot. One complete circulation of an electron's motion 
along the Aharonov-Bohm ring picks up the AB phase <fi 
generated by the enclosed magnetic flux, $. Note that 
Z r p {e;-4>) = Z;(e;cf>) and T p (e;0) ^_r p (e; -0)_for p = 

L, R, but r L (e; -<f>) + T R (e; -<P) = T L (e; <P) + T R (e; 0). 
These relations under the inversion of the magnetic flux 
leads to the (broken) AB phase symmetry in (out of) 
equilibrium, respectively. 

In the wide conduction band limit, all the rele- 
vant Green's functions of two electrodes become energy- 
independent except for the thermal functions. The trans- 
mission probability spectral function reduces to (see Ap- 
pendix for details) 



T(e) 



4", 



(1 + 7) 2 
1-7 



S^£-ImG d 

r L + r R 



= =- lmG r d {r L z r R 

i + ~/r L + r R dX R 



r^£}.(8) 



The first term is the transmission amplitude due to the 
direct tunneling. The second term is the expected form 
for transport through a quantum dot with the renormal- 
ized hopping rates Tl and T^. The third term is the 
interference between two current paths. We can further 
simplify the transmission probability spectral function 
T(e) as 



T(e) 



To 



To + 2lV.gTo(l - T ) cos ReG^ 

+f [T -,g(l -T cos 2 0)] lmG r d , (9a 

4 7 _ 4r L r R 



(1+7) 2 



(9b) 



Here T is equal to T L + T R = (T L + T R )/(1 + 7). Note 
that To is the transmission probability through the direct 
tunneling, and g is the maximum dimensionless linear 
conductance through a quantum dot in the absence of the 
direct tunneling. This dimensionless conductance also 
provides measures of the asymmetry in the coupling of a 
quantum dot to the left and right electrodes. 

In the linear response regime, we can expand the elec- 
tric and heat currents up to the linear terms of SV = 
Vl — V R and ST — Tl — T R . The transport coefficients 
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TABLE I: New parameters are tabulated in this table. 



(a) Zl(e; 0) = D ^ffi (e)| [T L + ^ R [g R {e)f + (z + z')g R {e)} 

(b) Z R (e;4>) = nr$0m [rR+ir L [gl(e)] 2 + (z + z*)gi(e)] 

(c) T L (e; 4>) = j§f^ [T L + yT R \g R (e)\ 2 + zg R {e) + z* g R (e)] 

(d) T R (e;4>) = [rR + jr L \gi(e)\ 2 +z*gi(e) + zgl(e)} 



Lij are defined by the relations 
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II 
Ql 



L21 L22 



Vl-Vr 
Tl-T r 



(10) 



and can be expressed in terms of the transport integral 

df 



I n (T) = IJde e"T(e) 



de 



(11) 



as L n = e 2 I , L 2 i = L 12 T = -eh and L 22 = h/T. The 
linear response conductance, G = limv^o dl/dV = Ln, 
is given by the equation, 



G = e 2 I (T), 



(12) 



which measures particle-hole symmetrical part of the 
transmission probability T(e) with respect to the Fermi 
level. The thermopower of a quantum dot in a 
two-terminal configuration can be found in an open 
circuit (I = 0) by measuring the induced voltage drop 
across a quantum dot when the temperature difference 
between two leads is applied. The thermopower is de- 
fined by the relation 



Vl-Vr 
hm — 



T, 



R 



1=0 



and can be expressed as 



S 



L\2 



h 



e ksTIo 



(13) 



(14) 



Here the constant fcs/e is approximately 86.17 /iV/K. 
In most metals, the thermopower is of the order of a few 
/iV/K. Since the integral I\ measures the first moment of 
energy e in the transmission probability, the thermopower 
probes the particle-hole asymmetric part in T(e) with 
respect to the Fermi energy. The thermal conductance k 
through an AB ring can be expressed in terms of these 
transport integrals as 



1 

T 



l2 ~T Q 



(15) 



The thermal conductance probes the particle-hole sym- 
metric part of the transmission probability. 



III. 



NON-CROSSING APPROXIMATION 



To study the thermoelectric effects of our model sys- 
tem, we have to find the Green's function, Gd, of a 



quantum dot. Due to the strong repulsive Coulomb in- 
teraction in a dot, the calculation of Gd is nontrivial 
and needs a many-body technique. For this purpose, 
we adopt the noncrossing approximation (NCA), which 
is a self-consistent diagrammatic method (for details, 
see Refsi22i2ii22i2ii2i). In our model system, the direct 
tunneling path renormalizes the Anderson hybridization 
functions and leads to a flux-dependent effective contin- 
uum band for the Anderson impurity (quantum dot), and 
the NCA integral equations are accordingly modified as 
shown below. 

The NCA has been very successful in studying the An- 
derson impurity models in the Kondo regime both in 
normal metals22*2i and out of equilibriumSSiSiSi. The 
zero-temperature analysis of the NCA self-energies re- 
veals that the NCA gives rise to the nonanalytic behavior 
of the Green's functions near the Fermi energy. However 
this nonanalytic behavior shows up in the finite temper- 
ature NCA below the pathological temperature T p . The 
value of T p can be estimated analytically 20 and T p <C Tk , 
where Tk is the Kondo temperature, so that the interest- 
ing Kondo effects in physical quantities can be computed 
down to temperatures far below Tk but still larger than 
T p . 

Due to the direct tunneling between the two electrodes, 
the effective continuum band for the Anderson impurity 
(quantum dot) is modified from two simple leads. In 
general, the self-energy for an electron on an Anderson 
impurity embedded in an effective electron continuum 
band can be decomposed into two terms, the one-body 
and many-body contributions, 



E d (e) = E c (e) + Etf(e), 



C + iS' 



(16a) 
(16b) 



Here Ej/(e) is the self-energy due to the on-site Coulomb 
interaction, while the first term E c comes from hopping 
into the continuum band. In our model system, the one- 
body contribution E c (e) is given by the Ea. l|A7b|l and 
can be written as 



Vd P G pP o(kt,k't')V pd . (17) 



p=L,R jjjf; 
p=LM ££/ 



The two auxiliary Green's functions, G p and G pp o, are 
introduced in Appendix ^ This one-body self-energy 
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(a) 




+ *V'\/\/V^\/\/\/\jb 




FIG. 2: NCA self-energy diagrams for (a) a pseudofermion 
and (b) a slave boson. The dashed (wavy) lines represent 
the propagators of pseudofermion (slave boson), respectively. 
The solid lines (look at Fig. llUl for definitions) are the effective 
continuum to the Anderson impurity model and represent all 
the multiple tunnelings between two leads (p = L,R). Open 
(solid) circle denotes the tunneling matrix between two leads 
(between a quantum dot and leads), respectively. 



shows up as the continuum band propagator in the Feyn- 
man diagrams of the pseudoparticle Green's functions, 
which will be defined below. 

In this paper, we take the limit of U — > oo and re- 
move the double occupancy in the highest-lying electron's 
state. Since we are interested in the Aharonov-Bohm ef- 
fect and the Kondo effect in the thermoelectric transport 
coefficients, the essential physics can be obtained using 
the infinite U NCA approach. In general the vertex cor- 
rection in the self-energy is needed to get the correct 
Kondo energy scaled when all three Fock spaces of an 
electron in a dot (empty, singly occupied, and doubly 
occupied configurations) are included. 

In the infinite U NCA, pseudoparticle Green's func- 
tions are introduced for the empty and singly occupied 
states, and the self-energy equations for these propaga- 
tors are solved self-consistently to second order in VdL 
and VdR. The empty state is represented by a slave bo- 
son operator b and singly occupied states (doubly de- 
generate) are denoted by the pseudofermion operator f a . 
The highest lying electron's annihilation operator in a 
dot can then be represented as 



(18) 



The NCA self-energy diagrams for a pseudofermion and 
a slave boson are displayed in Fig.|2|for our model system 
and their equations are 

£/(*,*') = ihG b (t,t')E c (t,t'), (19a) 
£ 6 (t,i') = -ihN s G f (t,t')X c (t',t), (19b) 

respectively. Here N s = 2 accounts for the two spin di- 
rections or the degeneracy of pseudofermion state. G b 
and G / are the Green's functions of the slave boson and 
pseudofermion operators and are defined as 



ihG b {t,t') = {T c b{t)b\t')), 
iW f (t,t') = (T c f a (t)p a (t>)). 



(20) 
(21) 



Projecting onto the physical Hilbert spaced, the lesser 
and retarded self-energies can be shown to be given by 
the equations, 

£/(<0 = / ^Gl(e-QX>(0, (22a) 

E?(6) = -J ^G<(e-()Xf(Q, (22b) 

££(6) = N.J !^(e + C)E<(0. (22c) 

^(e) = -N.J ^(e + OSrtC). (22d) 

These four equations are our NCA integral equations to 
be solved. The lesser and greater self-energies of a quan- 
tum dot due to the hopping into the two leads are 

£<(e) = 2r L (e;-0)/L(e) + 2rfl(e;-0)/ fl (e),(23a) 
E>(e) = 2r L (e;-0)/ L (e) + 2r fl (e;-0)/K(e),(23b) 

where the Anderson hybridization functions are given by 
the equations, 



r«(e; 



N L (e) 

N R (e) 
\Dr(e)\* 



[r L + 7 r R \g R (e)\ 2 

-z*g R (e) + zg a R {e)] , (24a) 



[r« + 7 r L | 5 £( e )| 2 

-zg r L (e) + z*gUe)}. (24b) 



Two Fermi-Dirac functions Jl and /r describe the ther- 
mal distribution in the left and right electrodes, respec- 
tively and f p = 1 — f p . The Anderson hybridization 
functions, Tl and Tr, are renormalized due to the direct 
tunneling term. In a wide conduction band limit, the 
renormalized Anderson hybridization becomes indepen- 
dent of energy variable and simplifies to 



\Pl,r + 1^r,l T 2y/~fT L r R sin ( 
+ 7? 



.(25) 



This reduction of couplings between the quantum dot 
and the two leads results in the smaller Kondo temper- 
ature than in the absence of the direct tunneling. In 
addition, the effective Anderson hybridization functions 
are explicitly dependent on the Aharonov-Bohm phase 4>- 
Note that the AB phase dependence of the Anderson hy- 
bridization functions is different in the above self-energies 
and in the equations of table |1] (c) and (d) denned in the 
expression for the transmission probability. When a fi- 
nite source-drain bias voltage is applied, Jl is not equal 
to /r, and the lesser and the greater self-energies do 
not remain invariant under the inversion of the magnetic 
flux, (j) — > —<j). Out of equilibrium, the AB phase symme- 
try is broken in the pseudoparticle self-energies and the 
Green's function of a quantum dot so that the transmis- 
sion probabilities are not the same under <p — > — <f>. In 
equilibrium, fx, = Jr and the AB phase symmetry under 
6 — > —<b is recovered. 
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IV. 



RESULTS 



Although we use the formalism presented in Ap- 
pendix lAl which treats electrodes with a general density 
of states (DOS), in actual numerical work a Lorentzian 
DOS is adopted, 



D 



iD 



(26) 



Here D is the bandwidth of two leads and is used as 
the energy unit (D = 1). The formalism for the flat 
DOS in the leads or the wide conduction band limit is 
useful in analyzing our numerical results. Throughout 
our numerical works, we use the following set of model 
parameters for a quantum dot. The energy level of a dot 
is chosen to be £d = —0.5D. The values of Tl and Tr 
are adjusted to satisfy both the total linewidth, Tl + 
Tr = 0.14D, and the chosen value of g (the asymmetry 
coupling factor). 

Since the effective Anderson hybridization function or 
linewidth of the quantum dot gets smaller due to the 
direct tunneling between two leads, the Kondo tem- 
perature Tk is also suppressed. In a wide conduction 
band limit, the linewidth is T = r^(— <f)) + Tr(— <j>) = 
(Tl + Tr ; ) / ( 1 4- 7) , which is independent of the AB phase 
4>. Tk can be estimated by the equation, 



Tk = Dy/c exp 



(27) 



where the dimensionless exchange coupling c is defined 
by the relation c = 2r/ 7r|ed| - One notable point is that 
the Kondo temperature is very sensitive to the value of 
the direct tunneling probability Jo- 
in the NCA approach to the Anderson impurity, the 
one-body hopping term (VdL and Van) is used as the ex- 
pansion parameter. Since only a subset of the self-energy 
diagrams is included up to the infinite order, the NCA un- 
derestimates the one-body contribution (£ c ) to the self- 
energy of an impurity site. Accordingly, the Fermi liquid 
relation^ is not satisfied at T = OK and the Kondo res- 
onance peak is exaggerated especially for the orbitally 
nondegenerate S = h Anderson model. This may lead 
to the violation of the causality relation well below Tk, 
or to the negative spectral weight in other renormalized 
Green's functions 2 ^! (the transmission probability T(e) 
near e = in our case). We can remedy this unphysical 
situation by exploiting the Fermi liquid relation 2 ^ for the 
self-energy of a dot, 



ImS^(e = 0,T = 0) = T. 



(28) 



This relation follows from the fact that the dot's self- 
energy due to the Coulomb interaction U satisfies the 
following relation^, 



ImE[/(e,T) oc e 2 + [irk B T\ A 



(29) 



when ksT, \e\ <C D,U. However the self-energy Y>nca 
of a dot computed from the NCA always satisfies the 
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FIG. 3: Temperature dependence of a Green's function of a 
dot near the Fermi energy, (a) The imaginary part or the 
spectral function develops the Kondo resonance peak near 
the Fermi energy with lowering temperature, (b) The real 
part develops a dip and peak structure with cooling. The 
temperature T is varied as T/T K = 100, 20, 4, 0.8, 0.16, 3.2 x 
10" 2 , 6.3 x 10" 3 , 1.3 x 10" 3 . The last four temperature curves 
cannot be distinguished with the naked eye. 



following inequality, 
— ImE 



NCA 



(e = 0,T=0) < r. 



(30) 



That is, the imaginary part of the NCA self-energy is 
always less than the Fermi liquid value, T. In fact, the 
NCA does a relatively good job in producing the energy 
and temperature dependence of the self-energy near the 
Fermi energy, but underestimates the one-body contri- 
bution to the dot's self-energy as noted above. Note that 
the one-body contribution is temperature- independent, 
and its form is known once the shape of the Anderson 
hybridization is given. We can make corrections in the 
one-body contribution to the NCA self-energy of a dot 
such that the above Fermi liquid relation is qualitatively 
satisfied. With this schem o 24 ! 27 , we remedy the unphysi- 
cal behavior of the NCA results at low temperatures well 
below Tk in our work. 

Transmission probability spectral function, T(e). The 
transmission probability shows the crossover behavior as 
a function of temperature due to the Kondo correlation 
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in the Kondo dot. At high temperatures above Tk , elec- 
trons flow from the left lead to the right one predomi- 
nantly via the direct tunneling since the tunneling of an 
electron from the leads to the dot is prohibited by the 
strong Coulomb repulsion (Coulomb blockade). In this 
case, the AB effect and the Fano interference in the trans- 
mission probability are weak. As the system is cooled be- 
low Tk, the Kondo resonance peak in the quantum dot 
develops close to the Fermi level and provides a new chan- 
nel for the flow of electrons though the dot. The newly 
opened current path interferes with the direct tunneling 
path leading to Fano interference and a strong AB effect 
in the transmission probability. 

The general structure of the transmission spectral 
functions near the Fermi level (e = 0) can be read off 
from the equation l|$af). As shown in Fig. ImG^ de- 
velops the Kondo resonance peak with its width of the 
order of Tk near e = while ReG^ varies very rapidly 
over the energy scale of Tk near e = with a dip just 
below e = and a peak above e = 0. The overall shape 
of the transmission spectral function T(e) is determined 
by the value of the AB phase <f> and the sign of A c [see 
Eq. (gjjfl, 



A,. 



T -.9(1 -T cos 2 , 



(31) 



A typical Fano interference pattern - a dip and peak 
structure - is expected when cos0 ^0. At cos0 = 0, 
T(e) has a dip(peak) resonance structure if A c > (<)0, 
respectively. Furthermore, we have one exact relation at 
zero temperature for the transmission probability at the 
Fermi level (e — 0), T(0) — g, which derives from the 
Fermi liquid relation 2 ^ for the Green's function of a dot, 
ImGS(O) = -i/r. 

In equilibrium, T(e, (f>) = T(e, — 4>) or the transmission 
probability remains invariant under the inversion of the 
magnetic flux, $ — > — <f>. This AB phase symmetry or 
Onsager relation can be deduced from the expression of 
T(e) in Eq.JSJ or Eq. l9a|) . If G r d (e) remains invariant un- 
der the transformation <j) — > —<f>, so is the transmission 
probability. In equilibrium, the Fermi-Dirac thermal dis- 
tribution functions in two leads are identical or /l = Jr. 
G d (e,—(f>) — G d (e, since the total Anderson hybridiza- 
tion Ti(e; —cj>) + Tn(e; —cj>) is an even function of </> and 
is reduced to T = (Tl + Tr)/(1 +7) in a wide conduc- 
tion band limit which does not depend on the AB phase. 
Due to the AB phase symmetry and the periodicity in <f>, 
we have only to consider T(e) in the range, < (j> < n. 
Out of equilibrium or when a finite source-drain voltage 
is applied, the two thermal distribution functions are not 
identical so the AB phase symmetry is broken in Gd and 
T(e). 

When the direct tunneling is weak, so is the Fano in- 
terference, and the Kondo-related resonance peak per- 
sists for the entire range of the AB phase. As the di- 
rect tunneling is increased, the overall magnitude of the 
transmission probability is also increased, and the Fano 
interference becomes more pronounced. 
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FIG. 4: Temperature dependence of the transmission proba- 
bility near the Fermi energy for the asymmetrically coupled 
dot. The model parameters are chosen as the asymmetrical 
factor g = 0.5 and the direct tunneling probability To = 0.3. 
The AB phase is varied as (a) A = 0°, (b) <f> = 90°, (c) 
A = 180°. 



When g = 1 or the quantum dot is coupled symmetri- 
cally to the left and right leads (r^ = T^), the transmis- 
sion probability T(e) for the case of T = 0.5, is presented 
in the short paper—. We briefly summarize the shape 
of T(e) in this case. When cos</> = 0, the transmission 
probability has a Kondo-related resonance peak near the 
Fermi level. When cos cj> 7^ 0, the typical shape of a Fano 
antiresonance shows up in the transmission probability. 
When < cos0 < 1, a dip (peak) structure develops be- 
low (above) the Fermi level (Ep), respectively. On the 
other hand, a dip (peak) structure arises above (below) 
Ep, respectively, when — 1 < cos0 < 0. 

In Figs. 21 and [SI we present the computed T(e) when 
the couplings are asymmetrical or g < 1. The AB phase 
dependence of T(e) is displayed in Fig. 0] when g = 0.5 
and To = 0.3. Comparing the curves in Fig. 0] and in 
Fig. 2 in our short pape»i&, the overall shape of T(e) 
is not qualitatively different between asymmetrically and 
symmetrically coupled dots when cos0 7^ 0. When </> = 
[see Fig. 01(a)], T(e) shows a typical Fano interference 
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FIG. 5: Temperature dependence of the transmission prob- 
ability near the Fermi energy for the asymmetrically cou- 
pled dot when 4> ~ 90°. The asymmetrical coupling factor 
is g = 0.5 and the direct tunneling probability is varied as 
(a) T = 0.1, (b) To = 0.5, (c) T = 0.9 from top to bottom 
panels. 



pattern, a dip structure below e = and a peak above 
e = 0. With further increasing <f> for 7r/2 < <j> < it, the 
Fano resonance pattern is inverted with respect to e = 
compared to the case < <f> < tt/2. For a symmetri- 
cal Kondo resonance peak at e = or when the energy 
structure in a quantum dot is particle-hole symmetric, 
exact inversion is expected, T(e, <f) — T(— e, 7r — 0), for 
< <fi < 7r/2. In our case of the infinite U calculation, the 
doubly occupied configuration in the dot is removed from 
our consideration, and particle-hole symmetry is broken. 
Due to this broken particle-hole symmetry, the spectral 
function of the dot has more spectral weight in the elec- 
tron excitations (e > 0) than the hole excitations (e < 0) 
and the relation T(e, <j)) = T(-e, ir - (p) for < (j) < tt/2 
does not hold in our case, which can be seen by compar- 
ing the curves in Fig. 21 (a) and (c). 

When cos <f> = or <f> = 90° , A c is equal to To — 9 and 
ReG^ is absent in T(e). The shape of T(e) near e = is 
solely determined by the sign of A c . A c remains always 
negative in the case of a symmetrically coupled dot so 
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FIG. 6: Electric conductance and AB effects. Model param- 
eters: g — 1 (symmetrically coupled dot). Direct tunneling 
probability (a) T = 0.1, (b) T = 0.5, (c) T = 0.9 from top 
to bottom panels. 



that the Kondo resonance peak persists in the transmis- 
sion probability spectral function over the entire range 
of the direct tunneling amplitude To. With increasing 
To, the Kondo-related peak becomes smaller and in the 
end is buried by the magnitude of T as T — > 1. Note 
that T(e) is less than or equal to one for all energies in 
our one transport mode model. The transmission spec- 
tral function reaches its maximum possible value g — 1, a 
unitary Kondo resonance tunneling. On the other hand, 
A c can change its sign from minus to plus in the case 
of an asymmetrically coupled dot. T(e) retains its peak 
structure near e = until A c < [See Fig. [S] (a)]. The 
maximum value of T(e) is set by the value of g. When 
A c = or T = g [see Fig. (b)], T(e) = T is a con- 
stant. The interference completely washes out the tun- 
neling probability through the quantum dot. The dip 
structure develops replacing the peak near e = as soon 
as A c becomes positive with increasing direct tunneling 
amplitude, T > g [See Fig. 03(c)]. The dip structure 
becomes more pronounced with further increasing To. 

Electric conductance, G(T, <fi). The temperature and 
AB phase dependence of the electric conductance is dis- 
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FIG. 7: Lorenz number and AB effects, g — 1 and Tb = 0.5. 



played in Fig. |H| for a symmetrically coupled dot. The 
temperature dependence of the electric conductance is 
determined by the values of Tq and <fi. When To is small, 
the electric conductance is monotonically increasing with 
lowering temperature for all AB phases. When To is 
large, G(T, 4>) is increasing for some range of and de- 
creasing for others (see Fig- El - 

The direct tunneling probability is a temperature- 
independent constant, while the transmission through 
the Kondo dot is strongly temperature dependent and 
shows crossover behavior as a function of temperature 
due to the Kondo correlation. Since the AB effects derive 
from the interference between two paths, the amplitude 
of the AB oscillation depends on the magnitude of the 
transmission coefficients for the two arms of the AB in- 
terferometer. A stronger direct tunneling amplitude, To, 
leads to a larger amplitude of the AB oscillation. 

The amplitude of the AB oscillation also shows 
crossover behavior as a function of temperature. The 
amplitude is weak (strong) at high (low) temperatures 
above (below) Tk, respectively. This crossover behavior 
derives from the Fano interference and the Kondo corre- 
lation in the quantum dot. Since the quantum dot lies 
in the Coulomb blockade regime at high temperatures 
above Tr-, the Green's function of the dot is negligibly 
small near the Fermi level so that the amplitude of the 
AB oscillation is weak. On the other hand, the Kondo 
correlation with lowering temperature below Tr- opens a 
new transport channel through the quantum dot and the 
amplitude of the AB oscillation is enhanced. 

Thermal conductance and Lorenz ratio, k(T, <p)/T and 
L(T,<fi). The diagonal transport coefficients, the electric 
and thermal conductances, are influenced by the Fano 
interference and the AB flux in the same way. The varia- 
tion of k(T,(/))/T versus T or is very similar to that 
of electric conductance G(T,(f>). The Lorenz number 
L(T, </>) is defined by the ratio of two diagonal conduc- 
tances 



L(T, 



k(T,<P) 



(32) 



TG(T,4>Y 

The variation of L(T) is displayed in Fig. [7| as a func 
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FIG. 8: Thermopower and AB effects, (a) g = 1 and To = 0.1, 
(b) g = 0.5 and T = 0.3 



tion of temperature T and the AB phase 4>. Our results 
obey the Wiedeman-Franz law satisfying the relation 

liniT^o L{T, <p) = ^--jf-, This relation is independent of 
the AB flux. The AB oscillation in the Lorenz number 
is strong in the crossover region or near the Kondo tem- 
perature. 

Thermopower, S. Since thermopower measures the 
average of the transmission probability weighted by the 
electronic excitation energy with respect to the Fermi en- 
ergy, it is a very sensitive probe of the particle- hole asym- 
metry in the transmission probability T(e). As shown 
in Figs. 0] and the shape of T(e) can be modulated 
through the Fano interference by the AB flux and the di- 
rect tunneling amplitude Tq. The sign of the thermopower 
can probe this structure. The sign of the thermopower 
is one of the experimental tools which can determine the 
main carriers of charge and heat. When more spectral 
weight lies in the electron excitations than in the hole 
excitations, the main carriers of charge and heat are 
electron excitations, and the sign of S is negative [see 
Figs-E^a) and[S{b)]. In the opposite case, the main car- 
riers are hole excitations, and the sign of S is positive 
[see Fig. 0Jc) and^b)]. In the particle- hole symmet- 
ric case (not shown here), electron and hole excitations 
carry the same amount of electric and heat currents. The 
signs of electric (heat) current are the same (opposite) for 
electron and hole excitations. The net result is that the 
thermopower is zero. The particle-hole symmetric case 
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can be realized when the quantum dot is located in the 
center of the Coulomb blockade diamond. 

Since the main part of our computed thermopower S 
in the AB interferometer is already reported in our short 
paper—, we briefly summarize the main results. We find 
that the sign of S can be changed by the energy structure 
in a quantum doliiS^ or the gate voltage capacitively 
coupled to the dot and by the AB flnsi^. The magni- 
tude of S is anomalously enhanced and is of the order of 
fc_B/e(« 86.17 fiV/ K) within the crossover region or near 
the Kondo temperature. For comparison, S is of order 
fiV/K in most of normal metals. The same enhancement 
of S is well known in heavy fermion systems and in the 
bulk Kondo systema 20 i 21 i 30 . 

When the direct tunneling amplitude is weak or T <C 
1, the sign of S remains invariant under the inversion 
of AB flux at temperatures near or less than Tk [see 
Fig.|H{a)]. This result follows from the fact that the Fano 
interference is weak and the Kondo-related peak persists 
in T(e) over the entire range of AB phase <f>. The sign 
change in S with AB flux can be realized when the direct 
tunneling To is increased [see Fig. EJb) and also Fig. 3 in 
our short pape»i&]. 

The amplitude of an AB oscillation in the thermopower 
is strongest within the crossover temperature region 
which is similar to the Lorenz number. The amplitude is 
weak (strong) when the direct tunneling is small (large), 
respectively [see Figs.IJa) and (b)]. 

The sign of S can be changed by varying the direct 
tunneling amplitude To. When To = 0-1 m Fig.OJa), the 
transmission probability has more spectral weight in the 
electron excitations than in the hole excitations. In this 
case, the thermopower is negative. In Fig. E^b) or when 
T = g = 0.5, T(e) is more or less particle-hole symmet- 
ric so that S ps 0. When the direct tunneling is increased 
or T) = 0.9 [Fig. E{c)], holes are the main carriers since 
the transmission probability is featured with more spec- 
tral weight below the Fermi energy. In this case, the 
thermopower is positive. 
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FIG. 9: Feynman diagrams for the mixed Green's functions 
which appear in the expression of current. Double wavy lines 
represent the Green's function of a quantum dot. The solid 
line with one arrow means the conduction electron propagator 
in the lead p — L, R. The solid line with two arrows is the 
auxiliary Green's function denned in Fig. HUf a'l. 



fined as the ratio of electric and thermal conductances, is 
fixed at the Sommerfeld value near zero temperature, but 
also shows AB oscillations near the Kondo temperature. 

The effect of the AB flux is more dramatic in the ther- 
mopower than in diagonal transport coefficients. In ad- 
dition to the AB oscillations in magnitude, even the sign 
of thermopower can be modulated by the AB flux. In the 
case of weak direct tunneling amplitude, the sign change 
is not possible under the AB flux since the Fano infer- 
ence is weak. When the direct tunneling amplitude is 
increased, the AB flux can change the sign of the ther- 
mopower. 

Most of studies in quantum dots have focused on the 
IV curves. Charge confinement is relatively easy with 
the use of the gate voltages, but the heat confinement is 
a more difficult job possibly because of the easy trans- 
fer of heat via other excitations. With more advances 
in nanotechnology, thermoelectric transport coefficients 
will provide additional information about the quantum 
transport in quantum dot systems. 
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V. SUMMARY AND CONCLUSION 

In this work we studied the effects of Fano interfer- 
ence and Aharonov-Bohm phase on the diagonal and 
off-diagonal transport coefficients in an AB ring with 
an embedded quantum dot in the Kondo regime. The 
transport properties shows crossover behavior from a 
high-temperature regime with no Kondo correlation to 
a low-temperature regime with Kondo correlation. The 
crossover behavior manifests itself in the amplitude of 
an AB oscillations in the diagonal transport coefficients, 
the electric and thermal conductances. The amplitude is 
small at high temperatures above Tk but becomes en- 
hanced below Tk by the Kondo correlation in the quan- 
tum dot. In the case of thermopower, the amplitude of 
an AB oscillation is strongest within the crossover region 
near the Kondo temperature. The Lorenz number, de- 
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APPENDIX A: DERIVATION OF EQUATION© 

In this section, we derive the expression of electric 
and thermal currents for general shape of the density of 
states(DOS) for two leads. Since the derivation of ther- 
mal current is the same as electric current, we derive the 
expression of electric current below in detail and discuss 
briefly thermal current at the end. 

Using the nonequilibrium Green's function method^, 
the electric current can be expressed in terms of the lesser 
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(a) Gn = 
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p 
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where two mixed Green's functions are denned by the 
equations 



(b) Z c = 



MGppo = 



p 

+ • — » 



p p 

> — o — » 



P 



P P 
» — o > 



ihG dp (k;t,t') = (Td a (t)cK k Jt')), (A2a) 
ihG w {kt,k't') = {Tc^J£)cl, p Jt')). (A2b) 



= Vdn 



° = T PP 



FIG. 10: Auxiliary Green's functions and one-body self- 
energy of a dot. Explaining the bare interaction vertex, the 
open circle means the direct tunneling matrix and the solid 
circle denotes the hopping amplitude between a quantum dot 
and two leads. 



mixed Green's functions, 



^ = TVY, I delm T LR G< L (k',k;e) 



h V 

kk' a 

2e 1 



= de Im 

ka 



V Ld G< L {k;e) , (Al) 



The Feynman diagrams of these two mixed Green's func- 
tions are shown in Fig. The mixed Green's function 
Gd P , shown in Fig.^a), is given by the equation 



G dp (k;t,t') = -^=V / dti G d (t,h)V dp G p (kit u ^) 
VV Jc 



The subscript C in the integral symbol denotes the 
Keldysh contour. The auxiliary Green's function G p is 
defined in Fig. HUf a) and is determined by the following 
Dyson-like equation 



G p (kt; k't') = G p (k;t,t') 6 S # + ^ £ 



k±k 2 



dh / dt 2 G p (kt;k 1 ti)T pp G p (k 2 ;ti,t2)T pp G p (k';t2,t'). (A4) 
c Jc 



r 



G p is the Green's function of the conduction electrons in 
the lead p when Tlr = = V dp . 

ihG p (k;t,t') = (Tc pi . a (t)cK k Jt')) Q . (A5) 

G d is the Green's function of a quantum dot, 

ihG d (k;t,t') = (Td a (t)dl(t')), (A6) 



and its self-energy (E d ) consists of two parts, the one- 
body contribution(S c ) coming from tunneling into two 
reservoirs and the many-body contribution (£u) due to 
the on-site Coulomb interaction U. Fig. HOf b'l gives the 
one-body contribution to the self-energy of a quantum 
dot. 



Z d {t,t') = £ c (t,t , ) + M*>*'), 
E c (t,0 



(A7a) 



^ Vd P G P (M, k't')V pd + ^Y.^Y.J dtl VdpGpfa, k^TppGpik'; t u t')V pd . (A7b) 

pkk' pkk 1 k\ 



S c includes all the multiple tunneling processes between two leads. Feynman diagrams of the mixed Green's function 
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G p p is shown in Fig. Elb) and G p p is given by the equation 



G pP (kt,k't r ) = Gppo(kt,k't') + - I d< i / d< 2 \G p (kt,kih) Vpd + Gppoikt,^) V p , 



k\.k-. 



C JC 



Vjjd 



xG d (t u t 2 ) [v dp G p (k 2 t 2 , k't') + V dp G pp0 {k 2 t 2 ,k't'] 
G pp0 (kt,k't') = dt 1 Gp(kt,k 1 t 1 )Tp P Gp{k'-t 1 ,t'). 

ki 

I 



(A8a) 
(A8b) 



Another auxiliary Green's function G pp q, defined in the electric current can be written as 
Fig. ITUT c'). is introduced to facilitate the algebra. In- 
troducing the wave- vector-summed Green's functions, 



Gpp(t,t) — — ^ ' Gpp(kt, k t ). 



Gd P {t,t') — 



(A9a) 
(A9b) 



2e 

7T 



A^ de Im [T Lfl G< L (e)] 



2e 



+— A, / de Im [V^G^e)] . (A10) 



(A9c) Here A, = 2 is two possible spin directions of conduc- 
tion electrons in the reservoirs. The wave- vector-summed 
Green's functions are determined by the following equa- 

(A9d) tions, 

I 



G dL {t,t') 
G RL (t,t') 



G L (t,t') 



dh D(t,h) [V dL G L (h,t') + V dR G RL o(ti,t')] , 



G RL0 (t,t')+ / dt x / dt 2 {G R (t,ti) V Rd + G R Lo(t,ti) VLd} 
JC Jc 

xG d (t u t 2 ) {V dL G L {t 2 ,t') +V dR G RL0 (t 2 ,t')} 

G RL o(t,t')+ [ dtx {G R {t,ti)V Rd + G R Lo{t,ti)V Ld } G dL (t 2 ,t'), 
Jc 

G L {t,t')+ [ dh [ dt 2 Gt(t, t±)TL R G R (ti,t 2 )T R LGL{t 2 , t'), 
Jc Jc 



G RL0 (t,t') = [ dtx G R (t,ti)T RL G L (ti,t') = f dt\ G R (t,ti)T R i,GL,(t\,t'). 
Jc Jc 



(Alia) 



(Allb) 
(Allc) 
(Alld) 



r 



Glr is obtained from the equation (jAllbjl by interchang- of the left and right reservoirs (summed over the wave 
ing L <-» R. vector) can be written as 

The retarded, advanced, and lesser Green's functions 

I 



V ^ e - e t ± id 

k P k 



= d( 



'e-(±iS 



kN p g r p > a (e), 



(A12a) 

G<(e) = 2irN p (e) f p {e) = 2irN p N p (e) f p (e) = 2nN p F p (e), (A12b) 

(Al2c) 

Here N p (e) — N p ■ N p (e) is the DOS of two leads and is normalized such that N p (0) — 1. Reduced Green's function 
g p ' a is introduced to simplify our algebra below and N p (e) = — Img^e). The auxiliary Green's functions, given by the 



G p (e) = 
G<(e) = 



Gppo ( e ) 

Gppo( e ) 



D r '°(e) 
1 

T - 
± pp 

D r > a (e) 

T - 
± pp 

\Dr{e)\- 

T - 
± pp 

\D r (eW 



= TTN; 



P £)r,a( c )' 



[G<(e) + |T Lfi | 2 |G£( e )| 2 G<( e )] 



G r p a (e)G^ a (e) 



T "" ^N L N R gT{e)g r A a {e), 



D r ' a (e) 



[G<(e)G a p (e) + G r p (e)G<(e)} 
2tt 2 N l N r [g%(e)F p (e) + g r p (e)F p (e)] . 



New variables are introduced, 



D^(e) = l-\T LR \ 2 GT(e)GT(e) = I - 7 g r £ a (e)g r A a (e), 
7 = ir 2 N L N R \T LR \ 2 , T p = irN p \V dp \ 2 . 



The electric current [Eg. ljA10|l ] can be written as 

II 



^Jde ImG< L (e), 

Gf L (e) = 2T LR G< L (e) + 2V Ld G< L (e) 
= 2T LR G RL0 {e) 

+2G^{VdLG L + V dR G RL0 }T LR {G R V Rd + G RL0 VLd} 
+2G r d {V dL Gt + V dR G< w }{T LR G r R V Rd + (1 + T LR G r RL0 )V Ld } 
+2Gj{V dL G L + V dR G RL0 }{T LR G R V Rd + (1 + T LR G RL0 )V Ld }- 



Note the following reduction of equations 



V dR G RLa = 


^ [VdL + V dR G 


Tl R G RL0 = 


1 




V dR G RL0 = 


i 


V dL G< + 


\D r \ 2 


G R wVLd = 


1 




\D r \ 2 



G| 
D' 



V 



(IL- 



d R u R 1 RL^> L 



* R 1 RL^'L v Ld 



where the renormalized tunneling matrices are introduced to simplify the equations. 

V dp = Vdp + VdpGp' a T pp , 

V pd = Vpd + T p pGp a Vp d . 

The current spectral function in the left lead can be written as 

4~; 



Gf L (e) 



\D r \ 



+ 



{g a L {e)F R {e) +g R ( e )F L (e)} 

2 Gj f -— a 



Fp^f {y dLG a L T LR G^V Rd + \T LR \ 2 G a L G r R ■ V dL G^V Ld j 



+ LE^P It { V <il g l V Ld + V d R G R T RL GlV Ld j 



G d ■ G a L \V dL \ 
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We also find the current spectral function in the right lead, 

47 



Gf R (e 



\D r \ 2 



{g a R (e)F L (e) + g r L (e)F R (e)} 



+ ■ 



2 Gl 



D r \ 2 D 
2 G 



+ 



\D r \ 2 D 
2 



f { v cIrGrTrlGlV Ld + \T LR \ 2 G a R G r L ■ V dR G^V Rd j 
I {v^RGtFu + V r dL GfT LR G R V r Rd } 



J^n G d ' G R \V dR \ 2 . 



(A19) 



Since only the imaginary parts are relevant to the expression of electric current, the current spectral function in each 
lead can be further simplified. 

lmG< L (e) = T (e) {f R (e) - f L (e)} + 2G< ■ T L (e) 

+ Wf 1 ^ { r dLGfV r Ld [l \T LR \ 2 GlG R ]+V dR G<T RL [Gl G L ]V Ld }} , (A20a) 
ImG<j(e) - T (e) {f L (e) - f R (e)} + 2G< ■ T R (e) 

+ ^p Im {^irG^vUi - \T LR \ 2 G R G a L ] + V r dL G<T LR [G R - G R ]V Rd )} , (A20b) 

where we introduced the transmission probability T for the direct tunneling and the renormalized Anderson hy- 
bridization r„ 



Me) 



r P (e) = 



47 



-N L (e)N R (e), 



*N p {e)\V a a 
\ D r(eW 



To simplify the algebra, we write 



ImG<(e) = To(e){/ p -(e)-/ p (e)} + 2G<-r p (e) + 2 £ A pq G<{e), 

q=L, R 



All = 
A LR = 
A RL = 
A RR = 



\D 



r\2 



^{^V r dL V r Ld [l-\T LR ?G 



\D r 



Gr 



—I m \-±V dR T RL [Gl-Gl]V Ld 



D 

G d 



737p Im \ ^ V il t lr[Gr - G R ]V Rd 



\D r 



g: 



-^Im \ -^V dR V Rd [l \T LR \ z G R Gl] \ . 



(A21a) 
(A21b) 

(A22a) 

(A22b) 
(A22c) 
(A22d) 
(A22e) 



Using the current conservation in a steady state, II + I R — 0, the expression of G d can be obtained in terms of the 
retarded and advanced Green's functions. 

1 [(All + A RL )Gf + (A LR + A RR )G<] . (A23) 



G< = 



T L (e)+T R (e) 

Inserting the lesser Green's function of a quantum dot into Gf L , we find the current spectral function, 
lmGf L (e) = T (e) [f R ( e )-f L (e)} 

[(TrAll - T L A RL )G< + (T R A LR - T L A RR )G<] . 



After some algebra, we can readily derive the following identities, 



r/.'-Ui - r, A llL - -^lJ^D a TrJVfl(e) V dL V Ld V dR V Rd , 



D r 



T R A LR -T L A RR = -_L^mgiD a n N L (e) V dL V r Ld V r dR V Rd . 



(A24) 

(A25a) 
(A25b) 
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Using these results, the current spectral function is reduced to a simple form, 



lmG< L (e) = T(e) [f R (e)-f L (e)}, 
T(e) = T (e)-lm[G r d Z LR }. 

Newly introduced parameters are given by the equations, 



(A26a) 
(A26b) 



Zlr 
Zl(e) 
Z r R (e) 

r*(6) 



AZIZ R 






T L (e)+T R (e)' 






nN L (e)V r dL V r Ld 




D r (e)\D r \ 


D r {t 


01^1 




N 


K(€) 


D r (e)\D r \ 


D r ( 


e)\D<"\ 


*N L (e)\V£ L \ 2 
\D r \ 2 


N L (e) 
\D r \ 2 


[r L + 


nN R (e)\V d a R \ 2 
\D r \ 2 


N R (e) 
\D r \ 2 


[Tr + 



-[T L+1 T R [g R {e)f + {z + z*)g R (t)} 
-[r R + 1 T L [gl(e)} 2 + (z + z*)gl(e)] 
[?l +jr R \g R (e)\ 2 + zg R (e) + z*g R {e)} , 
[T R + 7 T L \gl(e)\ 2 + z*gl(e) + zgl(e)] . 



(A27a) 
(A27b) 
(A27c) 
(A27d) 
(A27e) 



Here z = it 2 N lN R V d LTL R V R d = Vi^l^rb 1 ^ where <j> is the Aharonov-Bohm phase due to the magnetic flux threading 
through the AB ring. 

Thermal current can be derived in the same way. The only difference between electric and thermal currents is the 
multiplying factor in the integrands of the Landauer-Biittiker formula. The charge of an electron, — e, in the electric 
current is replaced by the energy of an electron, e — /zl, in the thermal current leaving the left lead. The thermal 
current is given by the expression, 



Ql - IJde (e-/i L ) T(e) [f L (e) - f R (e)}. 



(A28) 



Here T(e) is the transmission probability which we derived in the above. 



APPENDIX B: WIDE CONDUCTION BAND 
LIMIT 

In a wide conduction band limit, all the Green's func- 
tions related to the leads become energy-independent 
constants. Especially, gp' a (e) = Th N p (e) = 1, and 
D r ' a (e) = 1+7. The expression of the current can be 
highly simplified in this case. All the relevant Green's 
functions of two leads are given by the equations, 



G;> a (e) 
G<(e) 

G r p a (e) 



u ppO 



Gppo 



^N p f p {e), 



I+7' 
2ttN, 



(1 + 7) 2 



■k 2 N l N. 



1 + 7 
2iir 2 N L N R 
(1+7) 2 



" [f P (e)+lfp(e) 
T 



L ppi 



(Bla) 
(Bib) 

(Blc) 
(Bid) 
(Blc) 



Tpp[f P (e)-f p (e)}. (Blf) 



The parameters introduced in the previous Appendix are 
also simplified, 



Zr 



Tr 



1 



[T L -^T R -i(z + z*)], (B2a) 



(1 + 7) 2 

[T fl -7rL -*(* + **)]> (B2b) 

1 



[Vl+^Yr-^z-z*)}, (B2c) 



(1 + 7) 2 

j r ^[r R + 7 r L + i(z-z*)}. (B2d) 



Finally, the parameters in the transmission spectral func- 
tion become 
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r(i + 7 )i 



[Tl - 7^ - i(« + z*)] [T R - jTl - i(z + z*)] 



T Q - 0(1 - T cos 2 4>) + 2^ v / 5 T (l-T )cos ( /> 



(B3) 



Here T = T/(l + 7) and T = T L + T R . After some algebra, we find lmGf L to be lmGf L {e) = T{e) [f R (e) - f L {e)] 
with 



T(e) = - 47 . „ + _ 4rLl _^ ImG^ - 4 1 — I _ Im { G r d [T L Z R + T R Z* L ] } . 

(i + t) 2 t l + t r 1 + 7 r L + r fl 1 L R Li s 



(B4) 



r 



The transmission spectral function T(e) can be written To and g are the dimensionless conductance of direct 
in another form, tunneling and a resonant tunneling through a quantum 

dot, respectively. We note that the expression of T(e) 
in the wide conduction band limit was derived in the 
literaturei 16 i 31 i 32 



T(e) = T + 2ry.gT (l - T ) cos0 ReG r d 

+f [T -T cos 2 </>)] ImG r d , (B5) 

47 4r L r fl 



T» 



(l + 7 )2' 



(r 2 



r R ) 2 ' 



(B6) 
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